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>Method8  developed  for  analysis  of  intermittent  motion,  flexible  system  dynamics, 
design  sensitivity,  and  differential-algebraic  equations  are  summarized.  New 
techniques  for  dynamics  of  multibody  flexible  systems  include  mixed  vibration  and 
static  correction  modes.  Dynamic  design  sensitivity  analysis  methods  were 
developed,  using  singular  value  decomposition.  Numerical  Integration  methods  for 
mixed  differential-algebraic  equations  developed  are  based  on  singular  value 
decomposition  and  hybrid  generalized  coordinate  partitioning-constraint 
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Intenlttent  Motion 


A  reformulation  of  the  logical  function  method  for  representing  Impact 
events  was  developed,  using  a  coefficient  of  restitution  parameter  to 
establish  a  relationship  between  damping  and  elastic  effects.  Both  the 
logical  function  and  pieced  Interval  momentum  balance  methods  for  representing 
Intermittent  motion  effects  were  coded  Into  the  DADS  simulation  language  and 
tested  on  a  variety  of  problems.  The  pieced  Interval  technique  demonstrated 
excellent  numerical  performance.  In  addition,  the  pieced  interval  method  was 
reformulated  to  deal  with  constraint  addition  and  deletion,  for  analysis  of 
certain  classes  of  mechanisms.  Continued  development  and  application  for 
armaments  and  vehicles  was  pursued  under  an  exploratory  development  project. 
Flexible  System  Dynsmlcs 

Flexibility  effects  within  planar  and  three  dimensional  dynamic  systems 
were  represented  using  a  modal  coordinate  technique  to  describe  deformation  of 
bodies,  relative  to  a  body  fixed  reference  frame  that  Is  used  to  locate  the 
body  In  a  global  reference  frame  (1,  2,  3,  4].  A  general  formulation  is 
employed,  to  allow  the  user  to  obtain  normal  vibration  modes  and  stiffness  and 
mass  matrices  from  sources  such  as  experiment,  finite  element  analysis,  and 
analytical  techniques.  With  this  approach,  a  broadly  applicable  formulation 
has  been  obtained. 

Recently  developed  theories  of  static  correction  mode  methods  In 
transient  structural  dynamics,  developed  primarily  for  aerospace  and 
automotive  applications,  were  adapted  for  inclusion  in  formulation  for 
geometrically  nonlinear  system  dynamics  with  flexible  bodies  [5,6].  Test 
calculations  for  dynamics  of  mechanisms  and  machines  with  flexible  components 
were  carried  out  using  the  theory,  formulation,  and  experimental  computer  code 
[7,  8,  9).  Excellent  numerical  results  were  obtained  using  static  correction 
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mode  theory  from  structural  mechanics  to  account  for  local  deformations  that 
are  Induced  by  large  reaction  forces  that  occur  In  kinematic  joints.  Examples 
Created  Indicate  that  superior  results  can  be  obtained  using  the  static 
correction  mode  theory  In  conjunction  with  the  more  conventional  normal  modes 
of  vibration  of  structural  components.  Two  of  Che  examples  showed  that 
computing  cost  for  comparable  accuracy  can  be  as  much  as  three  orders  of 
magnitude  lower  If  static  correction  modes  are  used  In  conjunction  with  normal 
modes  of  vibration,  as  opposed  to  using  exclusively  normal  modes  of 
vibration.  These  results  were  not  totally  anticipated,  but  Indicate  the  value 
of  the  approach  for  future  applications  to  dynamics  of  high  performance 
robots,  vehicles,  aircraft,  space  structures,  and  other  mechanical  systems  In 
which  forces  transmitted  by  kinematic  connection  between  flexible  bodies  are 
Important. 

Design  Sensitivity  Analysis 

The  adjoint  variable  method  for  design  sensitivity  analysis  of  dynamic 
systems  was  extended  to  the  generalized  coordinate  partitioning  formulation, 
for  both  first  and  second  order  design  sensitivity  analysis  [10,  11].  Test 
calculations  with  moderately  simple  mechanisms  showed  that  accurate 
derivatives  of  dynamic  response  measures  with  respect  to  design  parameters  can 
be  obtained  and  used  for  dynamic  system  optimization.  Large  amounts  of  data 
accumulated  during  forward  Integration  In  time,  however,  must  be  stored  for 
backward  Integration  of  the  adjoint  equations.  This  factor  lead  to 
investigation  of  alternate  methods  for  numerical  Integration  of  the  equations 
of  motion,  to  minimize  the  amount  of  data  that  must  be  stored  for 
implementation  of  the  algorithm. 

A  numerically  efficient  adjoint  variable  design  sensitivity  analysis 
formulation  was  developed,  using  singular  value  decomposition  for  numerical 
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Integration  for  both  state  and  adjoint  equations  [12,  13,  14].  Symbolic 
computation  was  employed  to  generate  the  many  derivatives  required  with 
respect  to  design  variables  that  are  required  for  design  sensitivity  analysis, 
allowing  generality  In  user  definition  of  paraiaeters  to  be  declared  as  design 
variables.  The  formulation  was  Implemented  with  singular  value  decomposition 
and  tested  on  moderate  sized  examples  and  demonstrated  excellent  accuracy 
[15]. 

Numerical  Integration  of  Differential-Algebraic  Equations 

Development  of  singular  value  decomposition  for  defining  an  optimized  set 
of  Independent  generalized  coordinates  has  demonstrated  excellent  numerical 
stability  for  Integration  of  the  equations  of  motion  and  adjoint  equations  for 
design  sensitivity  analysis  [12,  13,  14,  IS]*  In  addition  to  much  Improved 
numerical  performance,  clear  physical  insights  have  been  gained  Into  the 
nature  of  composite  generalized  coordinates  that  are  defined  by  singular  value 
decomposition.  It  Is  shown  that  composite  generalized  coordinates  span  the 
tangent  hyperplane  to  the  constraint  surface,  associated  with  kinematic 
constraints  of  the  system. 

A  hybrid  constraint  stabilization-generalized  coordinate  partitioning 
method  was  developed  that  demonstrates  superior  computational  efficiency  and 
stability,  as  compared  to  methods  used  earlier  to  solve  mixed  differential- 
algebraic  equations  [16].  Numerical  criterion  for  111-conditlonlng  of  the 
constraint  Jacobian  with  respect  to  dependent  generalized  coordinates  was 
developed  and  Implemented,  as  a  positive  means  for  controlling  generalized 
coordinate  selection.  Numerical  experimentation  with  challenging  test 
problems  demonstrated  that  the  new  hybrid  algorithm  successfully  Integrates 
mixed  differential-algebraic  equations  that  are  Intractable  by  prior  methods. 
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